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short time scales to finite values of deterministic chaos at longer time scales for data 
generated from a Langevin equation that has a strange attractor in the limit of zero 
noise. Our results suggest that the crossover occurs at such short time scales that 
there is little chance of finite-dimensional Brownian noise being incorrectly identified 
as deterministic chaos. 
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A significant recent effort in nonlinear dynamics has been devoted to analyzing time 
series so as to determine whether they are chaotic or stochastic ]]]. Chaotic time series 
arise from deterministic strange attractors whose underlying geometric structure in state 
space can be used to improve short term forecasting [|I|], to remove noise or to stabilize 
unstable periodic orbits by small perturbations of parameters 0. Stochastic time series 
such as Gaussian white noise or Brownian motion can not be obtained from any finite set 
of deterministic equations, and are less amenable to prediction and control. 

An early expectation was that fractal dimensions of time series would provide a tool 
for distinguishing chaotic from stochastic behavior. A small empirical fractal dimension 
would correspond to a deterministic chaotic time series while a large or diverging fractal 
dimension (as a function of embedding dimension) would correspond to a stochastic time 
series |3],|]]. This conjecture was supported by some numerical experiments, e.g., the work 
of Ben-Mizrachi et al || which showed a crossover between high-dimensional white noise 
and low-dimensional chaos for time series generated by a chaotic map with additive white 
noise. The crossover occurred as a function of distance between points in phase space, with 
large correlation dimensions v (increasing linearly with the embedding dimension) occurring 
at small distances and a finite fractal dimension (independent of embedding dimension) 
occurring at larger distances. 

Several recent papers |||7||| have pointed out that this conjecture — that a fractal di- 
mension can distinguish deterministic from stochastic mechanisms — is not strictly correct. 
Osborne and Provenzale |J showed empirically that stochastic time series with a power-law 
power spectrum, P{oj) oc u a , yield a finite fractal dimension D that is parameterized by 
the index a; the precise value of D depends on some details of the random phases used in 
generating the time series ||. Theiler has pointed out that such time series are not likely to 
be identified incorrectly as deterministic if observed in experimental data [|7j . This follows 
since time series with a power-law power spectrum have long correlation times that scale 
with the length of the time series, and so are easily identified in the Grassberger-Procaccia 
dimension algorithm || by rejecting points that are close in phase space and also close in 
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time ||T0 . 

In this Brief Report, we explore further the extent to which stochastic Brownian motions 
can be distinguished from chaotic behavior using dimension calculations. Our starting point 
is the observation that the calculation of Mizrachi et al does not include the important case 
of a noisy dynamical system that evolves continuously in time such as a Langevin equation. 
Unlike a noisy map, the paths generated by noisy vector Langevin equations have a finite 
fractal dimension in the limit of short-time scales, with Hausdorff dimension D = 2 |]6|,|TT 



Over longer but still finite time scales (long compared to the correlation time of the time 
series), low-amplitude noise simply blurs the structure of the strange attractor and a finite 
fractal dimension is found, close to the dimension of the unperturbed attractor. 

There is then the possibility of an interesting crossover from finite-dimensional stochastic 
behavior to finite-dimensional deterministic behavior for continuous-time noisy systems. In 
the following, we study this crossover numerically and find that it occurs at such short time 
scales (for fixed moderate amounts of noise) that it is unlikely that the two effects will be 
confused in empirical data. The results below also confirm some scaling predictions made 
by Theiler in his analysis of fractional Brownian motions [0]. 

In our numerical experiments, we focus our attention on a vector Langevin equation of 
the form: 

X'(t) = F(X) + eG(t), (1) 

where the vector field F(X) is smooth and where the vector G represents a delta-function 
correlated Gaussian white noise source, (Gi(t)Gj(t')) = SijS(t — f). In the limit of large 
noise strength e ^> 1, we can ignore the field F and we obtain a Brownian-noise process 



whose fractal dimension is D = 2 [11|. In the opposing limit of small noise strength e<l, 
the Langevin equation reduces to a deterministic set of ordinary differential equations. To 
study a possible crossover in fractal dimension as the noise strength e is varied, we choose 



the vector field F to be a four- dimensional chaotic model derived by Lorenz 

F 1 = X 2 (X 3 -X A )-X 1 + c, (2a) 



F 2 = X 3 (X 4 -X 1 )-X 2 + c, (2b) 
F 3 = X 4 (X X - X 2 ) - X 3 + c, (2c) 
F 4 = X 1 (X 2 -X 3 )-X 4 + c. (2d) 

For the parameter value c = 100, these equations have chaotic solutions with a fractal 



dimension D pa 3.3 [T2|]. This dimension is distinctly larger than the value for Brownian 
paths. 

We expect the Langevin equation Eq. to have at least four different dimension scaling 
regimes in the general case, with the size of each regimes depending on the noise strength e, 
the length of the time integration T, and details of the strange attractor. For zero noise 
strength e = 0, we expect a scaling regime with dimension D = 1 at short time scales and 
a regime of D = 3.3 at long time scales compared to the correlation time. The former 
corresponds to the locally linear geometry of any smooth path in phase space. When noise 
is present, a third scaling regime with D = 2 (for Gaussian noise) should be visible at suffi- 
ciently short time scales. This new regimes arises from the scale-invariant non-differentiable 
structure of a Brownian path. For fixed noise and for sufficiently long time series, a fourth 
scaling regime with D — > oo should appear, representing the fact that the path eventually 
fills all of space. 

We have integrated numerically the Langevin equation Eq. p] with the vector field Eq. 2 
for different noise strengths e and lengths of time series N. We used a simple second- 



order accurate stochastic Runge-Kutta algorithm |13], and have verified that our results 
are not sensitive to the choice of time step At. Estimates of the correlation dimension v 
were obtained by delay-embeddings of time series of the first variable X\{f) into spaces of 
successively higher dimensions e, with an embedding delay of order the empirical correlation 
t c pa 0.1. Although the state vector X is known directly from Eq. |I| so that embedding is 
not strictly necessary, scaling regimes are identified more easily since they are invariant 
under increasing embedding dimension e. In our calculations, we also make the simplifying 
assumption that there is no observational noise, a complication that should be examined at 



a future date. 

Our results are presented in Figures 1-3. For zero additive noise e = 0, Fig. [I] shows 
the local slope v (correlation dimension) of scaling curves obtained with an approximate 
Grassberger-Procaccia algorithm, in which a sample of 1000 points is chosen from the time 
series and then all distances are calculated between each of these points and every other 
point in the entire time series. The 1000 points were chosen to be approximately evenly 
spaced throughout the time series. This method yielded similar results to Theiler's box- 
assisted method [14]], in which all distances between points in phase space are calculated for 



points that are sufficiently close. We do not implement Theiler's idea |L0[ of throwing out 
state-space points that are close in time, since these points give valuable information about 
possible stochastic structure. The scaling curves are presented over three different ranges 
in the phase-space distance r since a single calculation could not span the various scaling 
regimes. For quite short time scales At = 10~ 9 , the correlation dimension v is effectively one, 
confirming the locally linear smooth phase-space path. For larger sampling rates At = 10~ 5 , 
one sees a transition from v = 1 to larger values, but there is no apparent scaling suggestive 
of the strange attractor. For still larger sampling rates, now comparable with the typical 
time scales of the attractor, one sees a convergence of the local slope for different embedding 
dimensions around the known value D pa 3.3. This scaling regime is rather small, spanning 
just one order of magnitude in r. 

Similar calculations at the noise strength e = 1.0 are given in Fig. |2]. There is a new 
scaling regime with v pa 2.0 at the smallest time scales in panel (a). This indicates the finite 
fractal dimension of Brownian motion. This scaling regime changes smoothly onto the one 
of Fig. |l|(a), of locally linear orbits with v = 1. At still larger time scales, we see that the 
correlation dimension no longer converges for different embedding dimensions e, which is 
the hallmark of a higher dimensional process. 

The most important point of these two figures is that the crossover from Brownian noise 
to deterministic behavior occurs at quite small time scales, and at quite small distances r 
between points in phase space. This conclusion is qualitative in that the size of the scaling 



regimes will generally not be known a priori; the sizes will depend nontrivially on the noise 
strength e and on the length of the time series N. 

The difficulty of identifying all the relevant scaling regimes in dimension plots from 
time series sampled with a fixed time step motivated us to explore alternative ways to see 
this same structure. Towards this goal, we used an algorithm given by Higuchi [T^Q] that 



estimates a length-based fractal dimension D of the graph x(t) of a curve. This dimension 
is more appropriate for stochastic, rather than deterministic, series with continuous curves 
giving a value D = 1, and self-similar stochastic curves giving a value D = 1.5, which is the 



value expected for Brownian motion ||16|| . For a chaotic time series and sampling rates much 
longer than the correlation time we expect and find that D = 2 since the time series acts as 
a plane-filling curve. (We will discuss this point further in a future publication.) 

For three different noise strengths e, Fig. ^| gives the local slope (Higuchi dimension D) 
of scaling relations as a function of time interval r; the symbol r is not quite the same as 
that used in the Grassberger-Procaccia plots and here directly means a time interval, rather 
than a space-state interval. In Fig. 0(a), we see a result that is equivalent to that already 
observed in Fig. 0(a), namely that there is a smooth functional dependence at small time 
scales, for which D = 1, and this crosses over quite clearly to a chaotic regime of D = 2 at 
time intervals of order the characteristic time of the attractor. 

At moderate noise strengths such as e = 0.01, there is a new scaling regime at the 
smallest time scales with D = 1.5, arising from the Brownian motion of the Langevin path. 
This crosses over to D = 1 at longer time scales and then to D = 2 at still longer time scales, 
indicating the different kinds of dynamic structure present. For still larger noise strengths 
such as e = 1.0, the scaling regime for Brownian noise increases at the expense of the linear 
regime. For extremely large noise strengths e > 1, the deterministic structure becomes a 
small perturbation of the Brownian noise and only the D — 1.5 scaling regime should be 
observable. 

From the results summarized in these figures, we suggest that continuous-time processes 
arising from nonlinear Langevin equations — a case that models dynamical experiments per- 



turbed by external noise — lead to easily separated fractal-dimension scaling regimes cor- 
responding to stochastic Brownian motion at short time scales and deterministic chaos at 
longer time scales. This generalizes the earlier calculations of Mizrachi et al and of Proven- 
zale and Osborne. 
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FIGURES 

FIG. 1. Plots of the Grassberger-Procaccia dimension v for a fixed zero noise strength e = 0, 
and for three sampling rates of time series: (a) At = 1CP 9 ; (b) At = 10~ 5 ; (c) At = 1CP 2 . The 
parameter r here indicates a distance in a e-dimensional embedding space. For each sampling rate, 
a series of length N = 10 5 was studied. The curves denoted by squares, crosses, and triangles 
correspond to embedding dimensions of e = 4, e = 7, and e = 10 respectively. 

FIG. 2. The same as for Fig. [l| except the noise strength has a finite value e = 1. There is 
a scaling regime in (a) around D = 2 which is typical of Brownian noise, while the deterministic 
scaling regime in (c) at larger time scales is no longer well-defined. 

FIG. 3. Plots of the Higuchi dimension D for numerical time series of the Langevin equation, 
Eq. [I], for noise strengths e of: (a) e = 0.0; (b) e = 0.01; and (c), e = 1.0. The parameter r indicates 
the time interval separating pairs of points in a time series . The curve in each panel is a union 
of data obtained from three time series of 10 6 points with sampling rates of At = 10~ 3 , 10 -6 , and 
10~ 9 , since a single time series was not adequate to span all the scaling regimes. 
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Figure [I] David A. Egolf and Henry S. Greenside, Phys.Rev.E, "Stochastic to Deterministic 
Crossover of Fractal Dimensions for a Langevin Equation" 




Figure § David A. Egolf and Henry S. Greenside, Phys.Rev.E, "Stochastic to Deterministic 
Crossover of Fractal Dimensions for a Langevin Equation" 




Figure ^ David A. Egolf and Henry S. Greenside, Phys.Rev.E, "Stochastic to Deterministic 
Crossover of Fractal Dimensions for a Langevin Equation" 



